马尾 松 HDR 基 因 克 隆 及 其 对 旱 与 盐 胁 迫 的 响应 
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(南京 林业 大 学 林木 遗传 与 生物 技术 省 部 共 建 教育 部 重点 实验 室 南方 现代 林业 协同 创新 中 心 ， 南 京 
210037) 
摘 要 : 干旱 和 土地 盐 渍 化 是 制约 林业 可 持续 发 展 的 重要 因素 ， 植 物 在 遭受 生物 或 非 生 物 胁迫 时 ， 会 在 叶 
片 释 放 凋 类 等 挥发 性 物质 。1- 凑 基 -2- 甲 基 -2-E- 丁 烯 基 -4- 焦 磷酸 还 原 酶 CHDRO Æ MEP 途径 的 末端 活性 酶 ， 
具有 提供 前 体 墓 类 物质 和 主要 限 速 作用 。 为 探究 马尾 松 HDR 基因 是 否 参与 干旱 和 盐 胁迫 条 件 下 的 胁迫 响 
应 ， 该 研究 克隆 了 马尾 松 HDR 基因 开放 阅读 框 ， 并 初步 分 析 了 其 生物 信息 、 组 织 特异 性 表达 水 平和 功能 。 
结果 表明 : (1) PmHDR 基因 编码 区 长 度 为 1 458 bp， 编 码 485 个 氨基 酸 ， 其 编码 蛋白 包含 LytB/IspH 基 
因 超 家 族 的 核心 序列 和 PLNO2821 多 功能 结构 域 ， 属 于 1- 羟 其 -2- 甲 基 -2-(E)- 丁 烯 基 -4- 焦 磷酸 还 原 酶 家 族 。 
(2) PmHDR 密码 子 使 用 偏好 性 较 弱 ， 偏 好 使 用 A/U 结尾 的 密码 子 ， 烟 草 、 拟 南 芥 与 酿酒 酵母 更 适合 作为 
其 异 源 表达 受 体 。 (3) qRT-PCR 结果 显示 ，PmHDR 基因 在 马尾 松 老 叶 中 表达 量 最 高 ， 其 次 为 幼 叶 、 幼 茎 
和 老 共 ， 在 根 中 表达 量 最 低 。 (4) 构建 基因 表达 载体 pBI121-PmHDR 并 转化 拟 南 芥 ， 转 基因 拟 南 芥 对 干 
旱 和 盐 胁 迫 表 现 出 更 强 的 抗 道 性 。 以 上 研究 结果 表明 PmHDR 参与 植物 干旱 和 盐 胁 人 迫 的 响应 和 调节 ， 该 结 
果 可 为 马尾 松 抗 逆 育种 提供 一 定 的 理论 支持 。 
关键 词 : 马尾 松 ，PmHDR， 帖 类 化 合 物 ， 基 因 克 隆 ， 胁 迫 响 应 
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Cloning of HDR gene in Pinus massoniana and its response to drought 


and salt stress 


WANG Jiawen, YAO Sheng, SU Huan, LIU Kexin, ZHU Peihuang, JI Kongshu* 
(Key Laboratory of Forest Genetics & Biotechnology of the Ministry of Education, Co-Innovation Center for 


Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing 210037, China) 


Abstract: Drought and land salinization are inhibiting factors for the sustainable development of forestry, and the 
plants will release volatile substances such as terpenoids during life or non-biological stress. In order to investigate 
whether HDR gene of Pinus massoniana is involved in stress response under drought and salt stress, the open 
reading frame of PmHDR gene was cloned, and bioinformation, tissue specific expression and preliminary function 
were analyzed. The results were as follows: (1) The coding region of PmHDR gene is 1 458 bp, encoding 485 
amino acids, and its encoded protein contains the core sequence of LytB/IspH gene superfamily and PLN02821 
multifunctional structural domain, which belongs to the l-hydroxy-2-methyl-2-(E)-butenyl-4-pyrophosphate 
reductase family. (2) The PnHDR codon use preference was weak, with a preference for codons ending in A/U. 
Nicotiana tabacum, Arabidopsis thaliana and Saccharomyces cerevisiae were more suitable as its heterologous 
expression receptors. (3) qRT-PCR showed that the PmHDR gene was most highly expressed in old needles of 
Pinus massoniana, followed by young needles, young stems and old stems and least expressed in roots. (4) The 
gene expression vector pBI121-PmHDR was constructed and transformed into Arabidopsis thaliana. The transgenic 
A. thaliana showed greater resistance to drought and salt stress. These results indicate that PmHDR is involved in 


plant response and regulation to drought and salt stress, and these results can provide some theoretical support for 
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Pinus massoniana stress-resistance breeding. 
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马尾 松 (Pinus massoniana) 为 松 科 松 属 高 大 乔木 ， 是 我 国 古 老 的 本 土 针叶树 种 ， 也 是 我 国 特 有 的 速生 
树种 之 一 ， 广 泛 分 布 于 秦岭 一 淮河 以 南 ， 云 贵 高 原 以 东 17 个 省 、 自 治 区 、 直 辖 市 ， 为 我 国 分 布 面积 最 广 
的 针 叶 林 树种 〈 李 爱民 等 ，2008) 。 随 着 全 球 气候 变 暖 ， 水 资源 短缺 问题 日 益 严重 ， 大 范围 持续 性 干旱 和 
土地 盐 渍 化 对 我 国 的 农林 生产 造成 了 严重 威胁 ， 使 植物 生长 受 限 ， 并 使 植物 体内 的 有 机 物 成 分 发 生 改 变 ， 
影响 其 价值 与 产量 (Pichler & Oberhuber, 2007; 张 金 凤 等 ，2021) ， 对 我 国 的 农林 生产 造成 了 严重 威胁 ， 
因此 提高 树木 对 干旱 和 盐 胁迫 的 抗 性 一 直 是 森林 培育 的 主要 目标 之 一 。 马 尾 松 作为 我 国 重 要 的 木材 、 采 脂 
树种 和 荒山 造林 先锋 树种 之 一 《〈 徐 向 华 和 丁 贵 术 ，2006) ， 具 有 较 强 的 抗 逆 性 ， 但 重度 干旱 胁迫 仍 显 
制 马尾 松 苗 高 和 生物 量 的 增长 ， 严 重 影响 马尾 松 的 生长 和 繁殖 〈 全 文选 和 丁 贵 杰 ，2017) ， 因 此 研究 马尾 
松 对 干旱 和 盐 胁迫 的 适应 机 制 至 关 重 要 。 
植物 在 遭受 生物 或 非 生物 胁迫 时 ， 会 在 叶片 释放 挥发 性 物质 ， 包 括 碳 烯 类 和 脂肪 酸 衍 生物 等 ， 这 些 物 
质 具 有 直接 和 间接 的 信号 调控 作用 CLaothawornkitkul etal., 2008) ， 能 够 及 时 激活 脱落 酸 CABAO . K 
HREP (MeJA) 、 油 菜 素 内 醋 BR) 等 植物 激素 的 合成 和 信和 号 转 导 通路 ， 调 节 相 关 基 因 表 达 ， 以 缓解 
或 适应 胁迫 环境 〈 阁 秀峰 等 ，2007; Obata & Fernie, 2012) 。 其 中 昔 类 化 合 物 响应 干旱 胁迫 的 程度 较 高 ， 
© 在 植物 抵抗 逆境 胁迫 的 过 程 中 起 到 重要 作用 (Rastogi etal., 2019) . HXW WEE UEBER 
~ 丰富 、 结 构 最 多 样 化 的 一 类 《〈 罗 永明 等 ，2003) ， 其 以 异 成 二 烯 为 基本 单位 ， 共 同 前 体 物 为 异 戊 烯 焦 磷 酸 
APP) 和 其 异 构 体 二 甲 烯 焦 磷 酸 (DMAPP) 。IPP 和 DMAPP 在 生物 体内 的 合成 主要 有 甲状 成 酸 途径 
(MVA pathway) 和 丙酮 酸 途 径 (MEP pathway) 两 种 途径 《〈 陈 晓 明 ，2018) 。 其 中 MVA 途径 主要 参与 倍 
FWAR, MEP 途径 主要 参与 单 帖 、 二 昔 和 类 胡 葛 卜 素 等 的 合成 (Rohmer, 1999) 。1- 羟 基 -2- 甲 基 -2- 
E- 丁 烯 基 -4- 焦 磷酸 还 原 酶 CHDRO Æ MEP 途径 的 末端 活性 酶 ， 催 化 1- 羟 基 -2- 甲 基 -2-E- 丁 烯 基 -4- 焦 磷酸 
(HMBPP) 发 生还 原 反 应 并 以 5:1 的 比例 生成 IPP 与 DMAPP 的 混合 物 ， 具 有 提供 前 体 莫 类 物质 和 主要 限 
= WEH EMS 2014) 。 
I. 早期 发 现在 大 肠 杆菌 中 ，HDR 基因 对 MEP 途径 有 限制 性 作用 (Cunningham et al., 20000 ; 在 番茄 
= (Lycopersicon esculentum) RERAMA (Arabidopsis thaliana) J] E MAARE ^E 2S8] h 
p< 素 的 过 程 中 ，HDR 转录 水 平 明显 增加 ， 说 明 HDR 表达 与 类 胡萝卜 素 的 积累 有 关 (Botella-Pavia et al., 
20100 ; 对 银杏 (Ginkgo biloba) 的 研究 显示 ， 转 入 HDR 基因 的 银杏 与 野生 型 银杏 相 比 ， 总 内 酯 含量 显著 
"e 提高 〈 张 雯 等 ，2008) ; HERR (Rauvolfia verticillata) 的 研究 发 现 ，MeJA、 乙 酰 水 杨 酸 CASA) 、 
© ABA、 紫 外 光 (UV) 等 激发 子 可 以 诱导 上 调 HDR 基因 的 转录 水 平 (Chen etal., 2010) ; ÆRE fN 
(Dendrobium officinale) 中 ，HDR 基因 的 表达 受 水 杨 酸 (SA) 和 ABA 的 调控 ， 但 不 受 MeJA 影响 〈( 吴 秋 
菊 等 ，2015) 。 这 些 已 有 的 研究 均 表 明 HDR 基因 在 MEP 途径 中 具有 重要 的 限制 调控 作用 ， 与 植物 类 胡 葛 
卜 素 、 内 醚 等 东 类 化 合 物 的 合成 密切 相关 ， 且 受 调控 非 生 物 胁迫 响应 的 ABA、SA 等 植物 激素 诱导 ， 说 明 
其 可 能 参与 植物 对 非 生 物 胁迫 的 响应 。 
在 马尾 松 中 ，MEP 途径 上 游 关 键 酶 1- 脱氧-D- 木 酮 糖 -5- 磷 酸 合成 酶 的 编码 基因 PmDXS 和 1- 脱氧 -D- 木 
酮 糖 -5- 磷 酸 还 原 异 构 酶 的 编码 基因 PmDXR 已 被 成 功 克 隆 ， 并 进行 了 功能 和 表达 分 析 。DXR 是 MEP 途径 
的 第 一 个 关键 酶 ， 对 PmDXS 基因 的 研究 结果 显示 ， 用 机 械 损 伤 、15% 渗 透 剂 、H20，、 乙 烯 利 、MeJA、SA 
等 6 种 非 生 物 胁迫 或 激素 处 理 马 尾 松 后 ，PmDX5S 基因 的 表达 量 在 短 时 间 内 均 有 所 上 调 ;， 同 时， 在 400 
mmol:L-! NaCl, 800 mmol- L! D- H Z E ( D-mannitol) . pH 5.0. pH 9.0 等 非 生 物 胁 迫 环境 下 ， 
TransB/pET28a-PmDXS 重组 菌 生长 状况 显著 优 于 对 照 ， 表 明 PmDXS 基因 在 干旱 和 盐 胁迫 中 具有 正 调控 效 
M (CER, 2021) 。 这 些 结果 均 说 明 PmDXS 基因 参与 非 生物 胁迫 应 答 ， 在 马尾 松 的 逆境 胁迫 响应 中 起 一 
定 作 用 。DXR 是 催化 MEP 途径 第 二 步 反 应 的 关键 限 速 酶 ， 对 PmDXR 基因 的 研究 结果 显示 ， 相 比 于 野生 
型 拟 南 草 ，PmDXR 转基因 拟 南 芥 的 生长 发 育 在 不 同 浓度 NaCl、SA、MeJA、D-Mannitol 的 胁迫 条 件 下 受 
到 的 抑制 较 小 〈 朱 灵芝 ，2021) ， 说 明 PmDXR RAEE- EEF E35; H FEBR NL. DUE, TE 


为 MEP 途径 的 末端 限 速 酶 HDR 的 编码 基因 ，PmHDR 极 有 可 能 参与 马尾 松 在 干旱 和 盐 等 非 生物 胁迫 的 响 
应 过 程 。 本 文 将 克隆 马尾 松 PmHDR 基因 序列 ， 对 其 进行 生物 信息 学 分 析 和 组 织 特异 性 表达 分 析 ， 并 构建 


NS 


阐明 马尾 松 丫 类 


基因 表达 载体 转化 拟 南 芥 ， 对 比分 析 转 基因 拟 南 芥 和 野生 型 在 干旱 和 盐 胁 迫 条 件 下 的 生长 情况 ， 拟 探讨 以 
下 问题 : CI) PmHDR 蛋白 理化 性 质 、 结 构 、 所 属 基因 家 族 及 物种 间 的 进化 关系 ; (2) PmHDR 基因 的 组 
织 特异 性 表达 模式 ; 〈3) PmHDR 是 否 参 与 干旱 和 盐 胁 迫 等 非 生 物 胁迫 环境 下 的 胁迫 响应 。 以 期 为 进一步 


化 合 物 形 成 过 程 关键 酶 基因 的 潜在 生物 学 功能 提供 帮助 ， 并 为 马尾 松 抗 逆 育种 提供 参考 。 


1 材料 与 方法 
1.1 材料 


2 年 生 马 尾 松 实 生 苗 和 拟 南 芥 种 子 ， 均 取 自 林木 遗传 与 生物 技术 教育 部 重点 实验 室 〈 南 京 林 业 大 学 ) 。 

大 肠 杆 菌 Trelief M 5a 感受 态 细胞 和 pClone007 Blunt 载体 购 自 北京 擎 科 生 物 科 技 有 限 公司 ; 根 癌 农 杆 
V GV3101 感受 态 细胞 购 自 上 海 唯 地 生物 技术 有 限 公 司 ; Plant Total RNA Isolation Kit, Gel DNA Extraction 
Mini Kit. Plant DNA Isolation, ChamQ SYBR Color qPCR Master Mix. ClonExpress II One Step Cloning Kit 
购 自 南京 诺 唯 赞 生物 科技 股份 有 限 公 司 ; BamHI Sac I... QuickCut Buffer JJ Ej Takara; PCR Master Mix. 


lst Strand cDNA 

For Agarose Gels 
试验 中 引物 

1.2 方法 

1.2.1 PmHDR 基 


Synthesis SuperMix JJ A 322€ ^E gH C ERO. 股份 有 限 公 司 ; DNA ladder, Loading buffer 
和 TAE Buffer 购 自 上 海 捷 瑞 生物 工程 有 限 公司 。 
合成 和 序列 测序 由 南京 擎 科 生 物 科技 有 限 公 司 和 上 海 杰 李 生 物 技 术 有 限 公司 完 成 。 


因 的 克隆 参照 南京 诺 唯 赞 FastPure® Plant Total RNA Isolation Kit 的 说 明 书 提取 2 年生 马 


喇 


松 总 RNA， 使 
测 RNA 的 完整 


SuperMix for qPCR 将 RNA 反 转 录 为 cDNA， 于 -20 °C 条 件 下 保存 备用 。 


] Nanodrop One 超 微量 紫外 分 光 光 度 仪 检测 RNA 的 浓度 和 纯度 ， 并 通过 琼脂 糖 凝 胶 电泳 检 
性 。 选 取 纯 度 高 且 完 整 性 较 好 的 RNA， 利 用 反 转 录 试 剂 盒 Ist Strand cDNA Synthesis 


根据 火炬 松 (Pinus taeda). 同 源 基因 序列 ， 使 用 DNAMAN 软件 设计 PmHDR 特异 性 引物 ， 分 别 为 


PmHDR-F 


C 5'-ATGCCTTCGAGCCTCAGCTTTGC-3' ) S PmHDR-R ( 5'- 


TACCGTCTGCAACGCCTCCTCAT-33) 。 以 道 转录 获得 的 马尾 松 cDNA 第 一 链 为 模板 进行 PCR 反应 ， 扩 
增 体系 总 体积 为 50 nL， 包 括 高 保 真 酶 20 AL， 上 游 和 下 游 引 物 各 1 pL，cDNA 2.5 pL，ddH20 25.5 uL. 
PCR 反应 程序 为 95 °C 预 变 性 3 min; 95 °C 变性 15 s，60 °CiEK 45s, 72 °C 延伸 1 min, 4 35 个 循环 ; 


72 °C 延伸 5 min 
回收 纯化 ， 实 施 


通过 菌落 PCR 筛选 阳性 克隆 送 至 南京 擎 科 生 物 科技 有 限 公司 测序 。 
1.2.2 PmHDR 基因 的 生物 信息 学 分 析 利用 在 线 软件 对 PmHDR 基因 及 其 


. PCR 扩 增 产物 经 琼脂 糖 凝 胶 电泳 检测 后 ， 使 用 Gel DNA Extraction 产物 回收 试剂 盒 进 行 
连接 克隆 载体 (pClone007) 并 转化 至 大 肠 杆 菌 Trelief TM So 感受 态 细 胞 中 ， 培 养 单 菌 落 ， 


Mie 


i f TR 白 进行 生物 信 REA T. RB 


体 网 址 见 表 1; 利用 DNAMAN 软件 对 PmHDR 和 蛋白 进行 同 源 性 比 对 ;运用 MEGA 7.0 的 邻接 法 构建 系统 发 


育 进化 树 。 
表 1 生物 信息 学 分 析 在 线 软 件 及 网 址 
Table 1 Websites of bioinformatics prediction 
在 线 软件 网 址 18 
Online software Website Applications 
BLAST https://blast.ncbi.nlm.nih.gov/Blast.cgi 序列 比 对 
Sequence alignment 
. 理化 性 质 分 析 
ExPASy-ProtParam https://web.expasy.org/protparam Physicochemical droperty 二 
. RIKEST 
ExPasy-ProtScale https://web.expasy.org/protscale Hydrophilic analysis 
Wy ERR 
NCBI-CD https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi 保守 结构 域 分 析 


Conservative domain analysis 


在 线 软件 


Online software 


Cell-PLoc2.0 


SOPMA 


rm 
Website 


https://npsa-pra 


http://www.csbio.sjtu.edu.cn/bioinf/Cell-PLoc-2 


bi.ibcp.fr/cgi-bin/ 


npsa automat.pl?page-npsa sopma.html 


Phyre2 


TMHMM 2.0 


SignalP 5.0 


NetPhos 3.1 


YinOYang 1.2 


NetNGlyc 1.0 


1.2.3 PnHDR 同 义 密码 子 偏好 性 分 析 利用 CodonW 1.4.2 软 伯 
A3s、U3s、C3s、G3s、 相 对 同 义 密码 子 使 月 


EMBOSS 中 的 CUSP 在 线 软件 计 香 
拟 南 芥 、 烟 草 (Nicotiana tabacum) ~ 


http://www.sbg.bio.ic.ac.uk/phyre2 


http://www.cbs.dtu.dk/service/.TMHMM 


http://www.cbs.dtu.dk/services/SignalP-5.0 


http://www.cbs.dtu.dk/services/NetPhos 


http://www.cbs.dtu.dk/services/YinOYang 


http://www.cbs.dtu.dk/services/NetNGlyc 


Hf (relative synonymous codon usage, RSCU) 、 
Ceffective number of codons, ENC) 、 密 人 码 子 适 应 系数 (codon adaptation index, CAI) 、 最 优 密码 子 使 
Z& (frequency of optimal codons, FOP) . GC 含量 和 密码 子 第 3 位 上 的 GC 含量 〈GC3s) 等 。 运 用 


密码 子 的 


cerevisiae) 的 基因 组 密码 子 使 
1.2.4 PnHDR 基因 的 表达 分 析 


cDNA 稀释 至 10 ng: uL. 


频率 数据 来 
前 取 2 年 生 马 


用 液 氮 研磨 ， 提 取 各 样品 总 RNA， 通 过 琼脂 糖 凝 胶 ! 
RNA 的 浓度 和 质量 ， 将 符合 要 求 的 RNA 反 转 录 成 czDNA。 使 月 


根据 克隆 所 得 PmHDR 基因 序列 ， 使 
PmHDR-F (5'-AGAAATCGCAGAGCAGAA-3') 和 q-PmHDR-R (5'-CAACGCCTCCTCATCCTT-3') ， 


TUA (GenBank: KM496535.1) 为 内 参 基因 ， 


大 肠 杆菌 


C Escherichia coli) BER i| 
密码 子 使 用 统计 数据 库 (http://www. kazusa.or.jp/codon) . 
尾 松 根 、 幼 叶 、 老 叶 、 幼 茎 、 老 茎 等 五 个 部 位 的 样品 ， 


pe 
Applications 


亚 细 胞 定位 预测 


Prediction of subcellular localization 


二 级 结构 预测 
Secondary structure prediction 
三 级 结构 预测 
Tertiary structure prediction 
跨 膜 结构 分 析 


Transmembrane structure analysis 


t S HK T 
Signal peptide analysis 
磷酸 化 位 点 分 析 
Phosphorylation site analysis 
O- 糖 基 化 位 点 分 析 
O-glycosylation site analysis 
N- 糖 基 化 位 点 分 析 
N-glycosylation site analysis 


EF 计算 PmHDR 密码 子 的 使 用 特性 参数 ， 包 


了 效 密码 子 


分 别 


s; 95 "C 变 性 Ss，60 °C 退火 延伸 30 s，40 次 循环 ， 分 别 设置 3 次 生物 学 重复 和 3 次 技术 重复 ，PCR 反应 


束 后 进行 熔 解 曲线 分 析 。 使 
差异 显著 性 。 


1.2.5 基因 表达 载体 构建 使 用 pBI121-GUS 质粒 作为 真 核 表达 载 


择 BamHI 和 Sac I Pj 6 


PCR 扩 增 ， 通 过 琼脂 糖 凝 胶 电泳 检测 产物 质 


段 。 


用 限制 性 核酸 内 切 酶 BamH I 和 Sac I X} pBII21-GUS 质粒 进行 双 酶 切 ， 将 酶 切 后 的 质粒 切 胶 


与 PnHDR 基因 片段 连接 ， 
挑 取 单 菌落 进行 PCR PB 


还 重组 载体 质粒 pBIL21-PmHDR s 


SPSS 26 软件 进行 数据 分 析 并 作 图 ， 采 月 


并 立即 转化 Trelief M So 大 肠 杆菌 ， 于 含 卡 那 霉 素 的 LB 
性 检测 ， 将 阳性 单 克隆 扩 繁 后 送 至 南京 擎 科 生 物 科技 有 限 公 司 进行 测序 ， 并 由 其 返 


邓肯 多 重 比 较 法 分 析 各 组 织 表 达 


) ”和 PmHDR-SacI-R 


LH 


收 后 


括 
" 


酵母 《Saccharomyces 


使 


B.C Nanodrop One 超 微量 紫外 分 光 光 度 仪 检测 各 
H Nanodrop One 测定 cDNA 浓度 ， 并 将 各 


] Primier 5.0 软件 设计 实时 荧光 定量 PCR (qPCR) 引物 q- 


以 


] SYBR Green 进行 实时 区 光 定 量 表达 分 析 。 反 应 体系 共 10 
uL. &15 cDNA 模板 1 nL， 上 下 游 引 物 各 0.2 uL, SYBR Green Real-time PCR Master Mix (2x)5 uL, 无 
水 3.6 uL. qPCR 反应 在 ABI StepOne Plus RICE Æ PCR 仪 上 进行 ， 反 应 程序 使 用 两 步 法 : 95 "C 预 变性 


菌 
30 


结 


=E. 
EEN 


本 ， 根 据 测 序 得 到 的 PmHDR 基因 序列 ， 选 
切 位 点 ， 设 计 去 除 终止 密码 子 且 含有 酶 切 位 点 序列 的 重组 引物 PmHDR-BamHI-F 
C 5'-ACGGGGGACTCTAGAGGATCCATGCCTTCGAGCCTCAGCTT-3' 
CGATCGGGGAAATTCGAGCTCCCGTCTGCAACGCCTCCT-3') 。 


( 5" 
以 马尾 松 cDNA. 第 一 链 为 模板 进 
量 和 大 小 ， 并 切 胶 回收 得 到 含 酶 切 位 点 的 PmHDR 基因 ORF 


行 
F 


, 


固体 培养 基 上 培养 过 夜 ， 


m 
mug 


1.2.6 遗传 转化 拟 南 芥 及 胁迫 处 理 用 重组 质粒 冻 融 法 转化 农 杆菌 ， 扩 大 培养 阳性 单 苗 落 ， 花 序 侵 染 法 转化 


拟 南 芥 。 收 取 拟 南 芥 种 子 ， 依 次 用 75% 乙 醇和 20% 次 氧 酸 钠 消毒 ， 均 匀 播 种 在 含有 50 mg'mL- Kan 抗生素 
的 MS 固体 培养 基 上 进行 抗 性 筛选 ， 将 其 中 长 势 良好 的 拟 南 芥 幼苗 移 栽 至 营养 土 中 继续 培 
生 型 和 转基因 拟 南 草 的 叶片 ， 液 氮 研 磨 ， 使 用 南京 诺 唯 赞 生 物 科 技 股份 有 限 公司 生产 的 FastPure@ Plant 


养 。 分 别 采 集 野 


DNA Isolation 提取 基因 组 DNA， 并 以 其 为 模板 ， 用 PmHDR-BamHI-F 和 PmHDR-SacI-R 做 引物 ， 进 行 


PCR 3$, jit PmHDR 转基因 拟 南 和 芥 。 重 复 盘 选 1 次 ， 采 集 拟 南 芥 种 子 。 


分 别 配制 含 50 mmol.LLNaCl 和 含 $ 10. 50. 100 mmol: L?! D-H FIHI MS 固体 培养 基 ， 将 野生 型 
和 转基因 拟 南 芥 消 毒 后 分 别 点 播 在 培养 基 上 ， 每 个 培养 基 播种 20 个 种 子 。 培 养 10 4 后， 观察 拟 南 芥 发 芽 
及 生长 状况 ， 从 每 个 培养 基 中 挑选 10 个 成 功 发 芽 的 拟 南 芥 ， 测 量 并 记录 其 胚 根 长 ， 使 用 RStudio 软件 对 根 


长 数据 进行 分 析 。 
2 结果 与 分 析 
2.1 PmHDR 基因 的 克隆 及 序列 比 对 


扩 增 获得 PmHDR 基因 中 间 片 段 序 列 ， 长 度 为 1 458 bp， 与 琼脂 糖 凝 胶 电 泳 结 果 相 符 〈 图 1) 。 
M 1 2 


1 500 bp —> 


500 bp —> 


M 代 表 DNA Marker，1 和 2 代表 PmHDR 基因 。 
M represents DNA Marker, 1 and 2 represent PmHDR gene. 
图 1 PmHDR 中 间 片 段 电泳 检测 结果 


Fig.1 Electrophoresis results of intermediate fragment of PnHDR 


通过 NCBI 的 Blastn 在 线 软件 ， 将 PmHDR 基因 序列 与 其 他 物种 进行 比 对 。 输 出 结果 显示 ， 共 有 5 种 


植物 基因 序列 与 PnHDR 相似 度 高 于 85%， 分 别 为 赤松 (Pinus densiflora) 、 火 炬 松 CP. taeda?) 、 白 云 杉 
(Picea glauca) 、 北 美 云 杉 CP. sitchensis) 和 日 本 落叶 松 (Larix kaempferi) ， 其 中 赤松 和 火炬 松 的 IDS1 


基因 与 PmHDR 基因 序列 相似 度 高 达 98% 以 上 。 


通过 NCBI 的 Blastp 在 线 软件 对 PmHDR 的 氨基 酸 序 列 进行 比 对 ， 结 果 显 示 ，PmHDR 与 赤松 、 火 炬 松 


的 氨基 酸 序列 相似 度 最 高 ， 达 到 97% 以 上 ， 与 北美 云 杉 、 日 本 落叶 松 、 银 杏 、 川 又 CMorus notabilis) ~ 
澳洲 坚果 (Macadamia integrifolia) ~ *H9t (Juglans regia) . PAZ (Ananas comosus) 等 物种 的 相似 度 也 
均 高 于 78%， 说 明 HDR 基因 在 进化 过 程 中 较为 保守 。 选 取 与 PmHDR 氨基 酸 序 列 相似 度 较 高 的 几 个 物种 ， 


通过 DNAMAN 软件 进行 氨基 酸 序列 比 对 〈 图 2) ， 发 现 PmHDR 的 氨基 酸 序列 与 其 他 物 


一 致 ， 均 含有 4 


个 保守 的 半 胱 氮 酸 残 基 活性 位 点 ， 这 些 活 性 位 点 可 能 参与 催化 过 程 中 铁 硫 键 的 形成 《Lu et al., 2008) 。 
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深蓝 色 区 域 代表 几 个 物种 的 氨基 酸 序列 完全 一 致 ， 粉 色 区 域 代表 毛 基 酸 序列 高 度 保守 ， 浅 蓝 色 与 白色 区 域 代表 毛 基 酸 序列 差 
异 较 大 。 全 代表 半 胱 氮 酸 残 基 保守 位 点 。 


The dark blue regions represent complete amino acid sequence identity for several species, the pink regions represent highly conserved 


amino acid sequences, and the light blue and white regions represent highly divergent amino acid sequences. A represents conserved 
cysteine residues. 
图 2 PnHDR 氨基 酸 序列 与 其 他 物种 的 同 源 性 比 对 
Fig.2 Amino acid homology alignment of HDR gene from Pinus massoniana and other species 
2.2 PnHDR 的 生物 信息 学 分 析 
22.1 蛋白 一 级 结构 及 理化 性 质 分 析 
PmHDR 编码 和 蛋白质 的 分 子 式 为 C2415H3s44Nes40744S19， 理 论 分 子 质量 为 54.55 kD， 理 论 等 电 点 为 5.98， 
由 485 个 氨基 酸 组 成 ， 共 包含 71 个 负电 荷 氨 基 酸 碱 基 和 65 个 正 电荷 氨基 酸 碱 基 。 该 蛋白 中 含量 最 多 的 氨 
基 酸 为 赖 氨 酸 ， 占 总 氨基 酸 的 9.3%; 含量 最 少 的 为 半 胱 氨 酸 ， 仅 占 1.4% (图 3: A) 。 和 蛋白 脂肪 系数 为 
80.99， 总 平均 亲 水 性 为 -0.447， 说 明 该 蛋白 是 亲 水 性 蛋白 (图 3: B); 不 稳定 系数 为 28.98， 说 明 该 蛋白 
为 稳定 蛋白 。 
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A. PmHDR 编码 蛋白 的 氨基 酸 组 成 ， B. PmHDR $i 388 EA KRKE 
A. Amino acid compositions of protein encoded by PmHDR; B. Hydrophilicity and hydrophobicity analysis of protein encoded by 
PmHDR. 


图 3 PmHDR 基因 编码 蛋白 质 的 理化 性 质 分 析 
Fig.3 Amino acid compositions of protein encoded by PmHDR 

2.2.2 保守 结构 域 的 预测 与 分 析 

PmHDR 的 氨基 酸 序 列 保守 结构 域 预测 结果 显示 ，PmHDR 基因 编码 蛋白 在 第 60 至 第 483 个 氨基 酸 残 
基 处 具有 PLN02821 多 功能 结构 域 ， 该 结构 域 属于 LytB/IspH 基因 超 家 族 (Accession: cl19123) ， 该 家 族 
成 员 可 将 1- 羟 基 -2- 甲 基 -2-(E)- 丁 烯 基 -4- 焦 磷酸 CHMBPPO 转化 为 异 戊 烯 焦 磷酸 OPP) 和 二 甲 烯 焦 磷酸 
(DMAPPO ， 由 此 可 预测 PmHDR 基因 编码 蛋白 质 具 有 还 原 HMBPP 的 功能 (图 4: AD 。 
N 2.2.3 亚 细 胞 定位 预测 
© 使 用 在 线 软件 Cell-PLoc2.0 和 IncLocator 对 PmHDR 的 亚 细 胞 定位 进行 预测 ， 结 果 均 显示 PmHDR 编码 
nu 蛋白 位 于 叶绿体 中 的 可 靠 性 较 高 。 
-一 2.2.4 蛋白 质 跨 膜 结构 及 信号 肽 分 析 
^ 利用 在 线 软件 TMHMM Server v.2.0 对 Pm HDR 编码 和 蛋白质 的 跨 膜 结构 及 跨 膜 方向 进行 预测 分 析 ， 结 果 
co ”显示 该 蛋白 不 存在 跨 膜 结构 ， 全 部 位 于 膜 外 (图 4: BO 。 
使 用 在 线 工 具 SignalP-5.0 Server 对 PmHDR 编码 蛋白 进行 信号 肽 预测 分 析 ， 结 果 显 示 其 不 含有 信和 号 肽 ， 
表明 其 为 非 分 泌 重 白 (图 4: CO. 
2.2.5 蛋白 质 磷酸 化 位 点 和 糖 基 化 位 点 分 析 
使 用 NetPhos 3.1 Server 在 线 软件 进行 磷酸 化 位 点 的 预测 分 析 ， 结 果 显 示 PmHDR 编码 蛋白 中 共有 47 
个 磷酸 化 位 点 ， 包 括 23 个 丝氨酸 CSO ，20 个 苏 氨 酸 T ，4 个 酷 氨 酸 CYO (图 4: DO 。 

分 别 使 用 YinOYang 1.2 和 NetNGlyc 1.0 Server 在 线 软 件 进行 O- 糖 基 化 和 N- 糖 基 化 进行 预测 分 析 ， 结 
果 显 示 PmHDR 编码 蛋白 中 共有 4 个 0- 糖 基 化 位 点 ， 分 别 为 第 4 (S) 34 (S) 、58 CDD 236 S) 个 
氨基 酸 (图 4: EO ;共有 3 个 N- 糖 基 化 位 点 ， 分 别 为 第 129、 第 172、 第 397 NAER (图 4: F) 。 
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A. PnHDR 基因 保守 结构 域 预测 ，B. PmHDR $138 
预测 ，E. 0- 糖 基 化 位 点 预测 ，F. N- 糖 基 化 位 点 预测 。 


序列 位 点 
Sequence position 


长 白 的 跨 膜 结构 预测 ，C. PnHDR 编码 蛋 


序列 位 点 


Sequence position 


的 信号 肽 预测 ，D. 磷酸 化 位 点 


A. Prediction of conservative protein sequences encoded by PmHDR; B. Prediction of protein transmembrane region encoded by 


PmHDR; C. Prediction of protein transmembrane region encoded by PmHDR; D. Prediction of phosphorylation sites; E. Prediction of 


O-glycosylation sites; F. Prediction of N-glycosylation sites. 


图 4 PmHDR 编码 蛋白 的 生物 信息 学 分 析 
Fig.4 Bioinformatics analysis of protein encoded by PmHDR 


22.6 蛋白 二 级 结构 分 析 
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PmHDR 基因 编码 蛋白 的 二 级 结构 预测 结果 显示 该 蛋白 中 含有 a 螺旋 211 个 ， 占 比 43.51%; 含有 无 规则 
卷曲 164 个 ， 占 比 33.81%; 含有 延伸 链 78 个 ， 占 比 16.08%; 含有 B 转 角 32 个 ， 占 比 6.60% (图 5) 。 


Aou b PD IH 4l 


图 5 Pn HDR 编码 蛋白 二 级 结构 预测 
Fig.5 Prediction of secondary structure of protein encoded by PmHDR gene 
22.7 蛋白 三 级 结构 预测 
采用 同 源 建 模 法 进行 PmHDR 编码 蛋白 的 三 级 结构 预测 (图 6) ， 结 果 显 示 PmHDR 编码 蛋白 与 同 源 模 
板 c3dnfB[1- 羟 基 -2- 甲 基 -2-(E)- 丁 烯 基 -4- 焦 磷酸 还 原 酶 ] 的 序列 相似 度 为 33%， 置 信 度 为 100%。 


图 6 PnHDR 编码 蛋白 三 级 结构 预测 
Fig.6 Prediction of tertiary structure of protein encoded by PmHDR gene 
2.2.8 系统 进化 树 的 构建 
通过 邻接 法 Cneighbor-joining) 将 PmHDR 氨基 酸 序 列 与 其 他 物种 进行 比 对 并 构建 进化 树 (图 7) ， 从 
7 可 以 看 出 ，PmHDR 与 裸子 植物 聚 在 一 文 ， 与 赤松 、 火 炬 松 同 源 性 更 高 。 
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图 7 PmHDR 系统 进化 树 


Fig.7 Phylogenetic tree of PrnHDR 


2.3 PnHDR 密码 子 偏 好 性 分 析 


2.3.1 GC、GC3s、ENC、CAI 和 FOP 分析 


使 用 CodonW 1.4.2 对 PmHDR 的 A3s, U3s, C3s, G3s, ENC, CAI, FOP 等 参数 进行 分 析 ， 结 果 分 


别 为 0.400、0.468、0.272、0.155、47.77、0.215 和 0.363。 其 中 ENC 常 月 
为 20~61，ENC 值 越 小 表示 基因 密码 子 偏好 怕 


其 值 范围 


好 性 较 弱 。CAI 为 编码 该 蛋 


为 最 优 密 码 子 和 其 同 义 密 码 子 的 比值 ， 二 者 值 介 于 0~1， 其 值 越 高 ， 表 示 密 码 子 人 
的 CAI 为 0.215，FOP 为 0.363， 均 更 趋 近 于 0，i 
H EMBOSS 中 的 CUSP 在 线 软件 计算 总 GC 含量 和 GC3s， 


平 越 高 。PmHDR 基因 


选择 偏好 性 弱 ， 表 达 水 平 较 低 。 此 外 ， 使 月 


用 偏好 性 显著 (Mensah et al., 2019) 。PmHDR 基因 
白 的 所 有 密码 子 相 对 于 这 条 基因 都 使 用 最 优 密码 子 的 情况 


ERIR, EUA ENC 值 小 于 35 表示 基因 
的 ENC 值 为 47.77， 说 明 PmHDR 基因 的 密码 子 选择 偏 
下 的 适应 系数 ，FOP 
i 好 性 越 强 ， 基 因 表 达 水 
步 说 明 PmHDR 基因 对 密码 子 的 


于 分 析 基 因 的 密码 子 使 用 偏好 ， 
密码 子 使 


结 


显示 PmHDR 基因 编码 氨基 酸 时 密码 子 GC 含量 和 GC3s 值 分 别 为 39.92% 和 34.77%， 说 明 密码 子 偏好 以 


A/U 结尾 。 
2.3.2 相对 同 义 密码 子 使 用 度 分 析 


在 PmHDR 基因 选择 的 密码 子 中 ， 有 28 个 密码 子 的 RSCU 值 大 于 1.0， 为 PmHDR 偏好 密码 子 〈 表 


2) , CCA, AGA, GGA, ACU 和 CUU 等 5 个 密码 子 RSCU 值 大 于 2.0， 偏 好 性 较 强 ， 其 中 编码 精 氨 酸 
CArg) 的 AGA 密码 子 RSCU (ÑX 2.70, fi P EU. m3 XN? (Pro) W CCC, CCG 和 编码 Arg 的 
CGG 密码 子 的 RSCU 值 为 0， 说 明 它 们 可 能 不 参与 PmHDR 基因 的 翻译 过 程 。PmHDR 基因 偏好 密码 子 
绝 大 多 数 G/C 结尾 的 密码 子 RSCU 值 较 低 ， 也 进一步 表明 PmHDR 基因 的 密码 子 偏好 以 A/U 结尾 ， 表 达 水 
平 较 低 。 


表 2 PmHDR 相对 同 义 密码 子 使 用 度 
Table 2 RSCU of PmHDR gene 
d ENS EBEN 
SRO SWT 个 数 ”频率 比例 | EAR amy 个 数 ”频率 比例 
Amino RSCU Amino RSCU : 
acid Codon Number Frequency Fraction acid Codon Number Frequency Fraction 
KAR GCA 1.87 14 28.807 0.467 RAR CCA 2.25 9 18.519 0.562 
Ala GCC 0.53 4 8.230 0.133 Pro CCC 0.00 0 0.000 0.000 
GCG 0.27 2 4.115 0.067 CCG 0.00 0 0.000 0.000 
GCU 1.33 10 20.576 0.333 CCU 1.75 7 14.403 0.438 
EMAR UGC 0.86 3 6.173 0.429 谷 氨 酰胺 CAA 0.93 7 14.403 0.467 
Cys UGU 1.14 4 8.230 0.571 Gln CAG 1.07 8 16.461 0.533 
RARR GAC 0.37 5 10.288 0.185 AR AGA 2.70 9 18.519 0.450 
Asp GAU 1.63 22 45.267 0.815 Arg AGG 1.80 6 12.346 0.300 
RAR GAA 1.27 28 57.613 0.636 CGA 0.30 1 2.058 0.050 
Glu GAG 0.73 16 32.922 0.364 CGC 0.60 2 4.115 0.100 
KNAAR UUC 0.44 4 8.230 0.222 CGG 0.00 0 0.000 0.000 
Phe UUU 1.56 14 28.807 0.778 CGU 0.60 2 4.115 0.100 
甘氨酸 GGA 2.06 16 32.922 0.516 丝氨酸 AGC 1.16 6 12.346 0.194 
Gly GGC 0.13 1 2.058 0.032 Ser AGU 1.74 9 18.519 0.290 
GGG 0.65 5 10.288 0.161 UCA 1.94 10 20.576 0.323 
GGU 1.16 9 18.519 0.290 UCC 0.39 2 4.115 0.065 
组 氨 酸 CAC 0.67 3 6.173 0.333 UCG 0.19 1 2.058 0.032 
His CAU 1.33 6 12.346 0.667 UCU 0.58 3 6.173 0.097 
FRAR AUA 0.75 8 16.461 0.250 HAR ACA 1.11 10 20.576 0.278 
Ile AUC 0.56 6 12.346 0.188 Thr ACC 0.22 2 4.115 0.056 
AUU 1.69 18 37.037 0.562 ACG 0.56 5 10.288 0.139 
HMAR AUG 1.00 12 24.691 1.000 ACU 2.11 19 39.095 0.528 
Met 
MARR AAA 1.02 23 47.325 0.511 CAR UGG 1.00 8 16.461 1.000 
Lys Trp 
AAG 0.98 22 45.267 0.489 AR GUA 1.03 9 18.519 0.257 
TAR CUA 0.34 2 4.115 0.057 Val GUC 0.23 2 4.115 0.057 
Leu CUC 0.51 3 6.173 0.086 GUG 1.26 11 22.634 0.314 
CUG 0.86 5 10.288 0.143 GUU 1.49 13 26.749 0.371 
CUU 2.06 12 24.691 0.343 KAR UAC 0.55 3 6.173 0.273 
UUA 0.34 2 4.115 0.057 Tyr UAU 1.45 8 16.461 0.727 
UUG 1.89 11 22.634 0.314 * UAA 0.00 0 0.000 0.000 
天 冬 酰胺 AAC 0.87 10 20.576 0.435 UAG 1.00 1 2.058 1.000 
Asn AAU 1.13 13 26.749 0.565 UGA 0.00 0 0.000 0.000 


Hk: 下 划 线 表示 该 密码 子 RSCU>1，* 表 示 终 止 密码 子 。 
Note: The underline indicates that the value of RSCU-1, * indicates the stop codon. 
2.3.3 Pn HDR 外 源 宿主 的 选择 

将 PmHDR 基因 的 密码 子 使 用 频率 分 别 与 拟 南 芥 、 烟 草 、 酿 酒 酵母 和 大 肠 杆 菌 等 模式 物种 的 基因 组 进 
行 比较 分 析 〈 表 3) ， 将 比值 大 于 等 于 2.0 和 小 于 等 于 0.5 作为 筛选 密码 子 偏好 性 差异 较 大 的 标准 ， 以 此 统 
计 差 异 个 数 。 与 2 种 模式 植物 的 基因 组 相 比 ，PmHDR 基因 的 绝 大 部 分 密码 子 使 用 偏好 性 差异 较 小 ， 与 拟 
南齐 和 烟草 基因 组 密码 子 的 使 用 频率 差异 较 大 的 分 别 有 16 和 14 个， 说 明 在 进行 PmHDR 遗传 转化 试验 时 ， 
烟草 和 拟 南 芥 均 可 作为 异 源 表达 的 受 体 。 与 2 种 微生物 的 基因 组 相 比 ，PmHDR 与 大 肠 杆 菌 密码 子 偏好 性 
差异 很 大 ， 使 用 频率 差异 较 大 的 密码 子 有 24 个 ， 而 与 酿酒 酵母 相 比 使 用 频率 差异 较 大 的 密码 子 只 有 14 个 ， 
密码 子 偏好 性 差异 显著 小 于 大 肠 杆菌 ， 说 明 酿 酒 酵母 更 适合 作为 PmHDR 微生物 内 异 源 表达 试验 的 受 体 。 
表 3 PmHDR 与 部 分 模式 物种 基因 组 密码 子 使 用 偏好 性 比较 


Table 3 Comparison of codon usage preference between PmHDR and other model organisms 


y 


i 


s 


HT 烟草 酿酒 酵母 大 肠 杆菌 


SERE 密码 子 i A H 基因 组 基因 组 * NH PmHDR PmHDR  PmHDR  PmHDR 
Amino PmHDR Arabidopsis Nicotiana Saccharomyces Escherichia 
2 Codon . di . /At /Nt /Sc /Ec 
acid thaliana tabacum cerevisiae coli 
(At) (Nt) CSc) (Ec) 
KNAAR UUU 28.81 21.80 25.10 26.10 22.20 1.32 1.15 1.10 1.30 
Phe UUC 8.23 20.70 18.00 18.40 15.90 0.40 0.46 0.45 0.52 
TAR UUA 4.12 12.70 13.40 26.20 13.80 0.32 0.31 0.16 0.30 
Leu UUG 22.63 20.90 22.30 27.20 13.00 1.08 1.01 0.83 1.74 
CUU 24.69 24.10 24.00 12.30 11.40 1.02 1.03 2.01 2.17 
CUC 6.17 16.10 12.30 5.40 10.50 0.38 0.50 1.14 0.59 
CUA 4.12 9.90 9.40 13.40 3.90 0.42 0.44 0.31 1.06 
CUG 10.29 9.80 10.20 10.50 51.10 1.05 1.01 0.98 0.20 
HZA AUU 37.04 21.50 27.80 30.10 29.70 1.72 1.33 1.23 1.25 
Ile AUC 12.35 18.50 13.90 17.20 23.90 0.67 0.89 0.72 0.52 
AUA 16.46 12.60 14.00 17.80 5.50 1.31 1.18 0.92 2.99 
MAR AUG 24.50 25.00 20.90 27.20 1.01 0.99 1.18 0.91 
Met 24.69 
谷 氮 酰胺 CAA 14.40 19.40 20.70 27.30 14.70 0.74 0.70 0.53 0.98 
Gln CAG 16.46 15.20 15.00 12.10 29.40 1.08 1.10 1.36 0.56 
EAR UGG 16.46 12.50 12.20 10.40 15.30 1.32 1.35 1.58 1.08 
Trp : 
Eu GUU 26.75 27.20 26.80 22.10 18.10 0.98 1.00 1.21 1.48 
Val GUC 4.12 12.80 11.10 11.80 14.80 0.32 0.37 0.35 0.28 
GUA 18.52 9.90 11.40 11.80 10.90 1.87 1.62 1.57 1.70 
GUG 22.63 17.40 16.70 10.80 26.20 1.30 1.36 2.10 0.86 
丝氨酸 AGU 18.52 14.00 13.30 14.20 9.40 1.32 1.39 1.30 1.97 
Ser AGC 12.35 11.30 10.00 9.80 16.00 1.09 1.23 1.26 0.77 
UCU 6.17 25.20 20.00 23.50 8.70 0.24 0.31 0.26 0.71 
UCC 4.12 11.20 10.20 14.20 8.90 0.37 0.40 0.29 0.46 
UCA 20.58 18.30 17.60 18.70 8.10 1.12 1.17 1.10 2.54 
UCG 2.06 9.30 5.30 8.60 8.80 0.22 0.39 0.24 0.23 
BAR CCU 14.40 18.70 18.70 13.50 7.20 0.77 0.77 1.07 2.00 
Pro CCC 0.00 5.30 6.60 6.80 5.60 0.00 0.00 0.00 0.00 
CCA 18.52 16.10 19.80 18.30 8.40 1.15 0.94 1.01 2.20 
CCG 0.00 8.60 5.00 5.30 22.40 0.00 0.00 0.00 0.00 
HAR ACU 39.10 17.50 20.30 20.30 9.10 2.23 1.93 1.93 4.30 
Thr ACC 4.12 10.30 9.70 12.70 22.80 0.40 0.42 0.32 0.18 
ACA 20.58 15.70 17.40 17.80 8.10 1.31 1.18 1.16 2.54 
ACG 10.29 7.70 4.50 8.00 15.00 1.34 2.29 1.29 0.69 
AAR GCU 20.58 28.30 31.20 21.20 15.40 0.73 0.66 0.97 1.34 
Ala GCC 8.23 10.30 12.50 12.60 25.20 0.80 0.66 0.65 0.33 
GCA 28.81 17.50 23.10 16.20 20.70 1.65 1.25 1.78 1.39 
GCG 4.12 9.00 5.80 6.20 32.30 0.46 0.71 0.66 0.13 
KAR UAU 16.46 14.60 17.80 18.80 16.50 1.13 0.92 0.88 1.00 
Tyr UAC 6.17 13.70 13.50 14.80 12.30 0.45 0.46 0.42 0.50 
HAR CAU 12.35 13.80 13.40 13.60 12.80 0.89 0.92 0.91 0.96 
His CAC 6.17 8.70 8.70 7.80 9.40 0.71 0.71 0.79 0.66 
RZ AAU 26.75 22.30 28.00 35.70 19.20 1.20 0.96 0.75 1.39 
Asn AAC 20.58 20.90 17.90 24.80 21.70 0.98 1.15 0.83 0.95 
赖 氨 酸 AAA 47.33 30.80 32.60 41.90 34.00 1.54 1.45 1.13 1.39 
Lys AAG 45.27 32.70 33.50 30.80 11.00 1.38 1.35 1.47 4.12 
RZA GAU 45.27 36.60 36.90 37.60 32.80 1.24 1.23 1.20 1.38 
Asp GAC 10.29 17.20 16.90 20.20 19.20 0.60 0.61 0.51 0.54 
RAR GAA 57.61 34.30 36.00 45.60 39.30 1.68 1.60 1.26 1.47 
Glu GAG 32.92 32.20 29.40 19.20 18.70 1.02 1.12 1.71 1.76 
半 胱 氨 酸 UGU 8.23 10.50 9.80 8.10 5.20 0.78 0.84 1.02 1.58 
Cys UGC 6.17 7.20 7.20 4.80 6.40 0.86 0.86 1.29 0.96 
甘氨酸 GGU 18.52 22.20 22.30 23.90 24.20 0.83 0.83 0.77 0.77 
Gly GGC 2.06 9.20 11.20 9.80 28.10 0.22 0.18 0.21 0.07 
GGA 32.92 24.20 23.20 10.90 8.90 1.36 1.42 3.02 3.70 
GGG 10.29 10.20 10.50 6.00 11.80 1.01 0.98 1.71 0.87 
HAR CGU 4.12 9.00 7.50 6.40 20.20 0.46 0.55 0.64 0.20 
Arg 


CGC 4.12 3.80 3.90 2.60 20.80 1.08 1.06 1.58 0.20 
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拟 南 芥 烟草 酿酒 酵母 大 肠 杆菌 


EN FE K] 2 Z Z 其 [K] 4] 
Aa 酸 密码 子 s x 基因 组 基因 组 B 组 .  PnHDR PmHDR PmHDR PmHDR 
Amino PmHDR Arabidopsis Nicotiana Saccharomyces Escherichia 
3 Codon . idi . /At /Nt /Sc /Ec 
acid thaliana tabacum cerevisiae coli 
CAt) Nt) CSc) (Œc) 
CGA 2.06 6.30 5.30 3.00 3.80 0.33 0.39 0.69 0.54 
CGG 0.00 4.90 3.70 1.70 6.20 0.00 0.00 0.00 0.00 
AGA 18.52 19.00 16.00 21.30 2.90 0.97 1.16 0.87 6.39 
AGG 12.35 11.00 12.20 9.20 1.80 1.12 1.01 1.34 6.86 
* UAA 0.00 0.90 1.10 1.10 2.00 0.00 0.00 0.00 0.00 
UAG 2.06 0.50 0.50 0.50 0.30 4.12 4.12 4.12 6.86 
UGA 0.00 1.20 1.00 0.70 1.10 0.00 0.00 0.00 0.00 


注 : 下 划 线 表示 密码 子 使 用 频率 差异 较 大 ，* 表 示 终 止 密码 子 。 

Note: The underline indicates that the codon usage frequency is quite different, and * indicates the stop codon. 

2.4 PnHDR 的 组 织 特异 性 表达 分 析 

通过 实时 荧光 定量 表达 技术 ， 分 析 PmHDR 基因 在 成 年 马尾 松 根 、 幼 叶 、 老 叶 、 幼 茎 、 老 茎 中 的 表达 
情况 。 结 果 表明 ，PmHDR 基因 在 马尾 松 根 中 表达 量 最 低 ， 在 老 叶 中 表达 量 最 高 ， 为 根 中 表达 量 的 122.2 fis, 
在 幼 叶 、 幼 蔡 、 老 茎 中 的 表达 量 分 别 为 根 的 21.2 倍 、6.5 倍 和 6.0 倍 《图 8) 。 
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不 同 小 写字 母 表 示 差 异 显著 (P<0.05) . 
Different lowercase letters indicate significant differences (P«0.05). 
图 8 PnHDR 基因 的 组 织 特异 性 表达 
Fig.8 Tissue specific expression of PnHDR gene 

2.5 胁迫 处 理 下 拟 南 芥 生 长 观测 

将 转基因 拟 南 芥 种 子 播种 在 含 Kana 抗生素 的 MS 固体 培养 基 上 进行 筛选 ， 并 选取 其 中 正常 生长 发 育 
的 植株 幼苗 转移 至 营养 土 中 继续 培养 2~3 周 后 进行 分 子 鉴 定 〈 图 9) 。 结 果 显 示 ， 被 选择 的 8 株 转基因 拟 
南 芥 中 有 6 株 检 测 结果 为 阳性 ， 因 此 采集 这 6 株 拟 南 芥 的 种 子 ， 开 展 后 续 实 验 。 


15000p Mg sooo —-— 5 weaf -—* 
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WT 为 野生 型 拟 南 芥 ，1-8 为 转基因 拟 南 芥 。 
WT is wild Arabidopsis thaliana; 1-8 is transgenic Arabidopsis thaliana. 
图 9 转基因 拟 南 芥 PCR 检测 
Fig. 9 PCR detection of transgenic Arabidopsis thaliana 

利用 D- 甘 露 醇和 NaCl 分 别 对 拟 南 芥 进行 干旱 胁迫 和 盐 胁 迫 处 理 ， 将 野生 型 和 PmHDR 转基因 拟 南 芥 
! 子 分 别 点 播 在 含有 不 同 浓度 D-H SEISURI NaCl 的 MS 固体 培养 基 上 培养 10 d， 观 察 并 记录 各 个 培养 基 
拟 南 芥 的 生长 发 育 情况 〈 图 10) 。 结 果 表 明 ， 在 正常 培养 条 件 〈 无 胁迫 处 理 ) 下 ， 转 基因 拟 南 芥 的 胚 根 长 
略 高 于 野生 型 拟 南 芥 ， 但 差异 不 显著 ; 与 之 相 比 ， 经 5S、10、50、100 mmol L^ D-H E EI 50 mmol L 
NaCl 处 理 的 拟 南 芥 生 长 受到 不 同 程度 的 抑制 ， 其 中 野生 型 拟 南 芥 平均 胚 根 长 分 别 降 低 35.49%、39.26%、 
53.59%、77.08%、77.06%， 转 基因 拟 南 齐 平均 胚 根 长 分 别 降低 11.03%, 7.40%, 21.8896. 59.0990. 20.1096. 
在 D- 甘 露 醇和 NaCl 胁迫 处 理 下 ， 转 基因 拟 南 芥 的 胚 根 长 显著 高 于 野生 型 拟 南 芥 ， 说 明 转 基因 拟 南 芥 受 D- 
露 醇和 NaCl 胁迫 影响 小 ， 对 干旱 和 盐 胁迫 表现 出 更 强 的 抗 逆 性 〈 图 11) 。 
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WT is wild Arabidopsis thaliana; PnHDR is transgenic Arabidopsis thaliana. The same below. 
图 10 不 同 处 理 拟 南 芥 生 长 状况 


Fig.10 Growth of Arabidopsis thaliana in different treatments 
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## 表 示 P<0.01，*** 表 示 P«0.001. 
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图 11 不 同 处 理 拟 南 芥 胚 根 长 


Fig.11 Root length of Arabidopsis thaliana radicle in different treatments 


3 讨论 与 结论 


苹 类 化 合 物 是 植物 次 生 代 谢 产 物 中 最 丰富 、 结 构 最 多 样 化 的 一 类 ， 在 植物 生长 发 育 和 抵抗 逆境 胁迫 的 
过 程 中 到 重要 作用 。1- 羟 基 -2- 甲 基 -2-E- 丁 烯 基 -4- 焦 磷酸 还 原 酶 〈(HDR ) 是 熙 类 化 合 物 合 成 中 MEP 途径 的 


末端 活性 酶 ， 上 共有 提供 前 体 向 类 物质 和 主要 限 速 作 用 。 为 探究 马尾 松 对 干旱 和 盐 胁迫 的 响应 ， 本 文 看 


EH 


EET PmHDR 基因 。 而 莫 类 化 合 物 生物 合成 MEP 途径 上 游 关 键 酶 的 编码 基因 PmDXS 和 PmDXR 已 被 成 功 

克隆 和 分 析 ， 发 现 二 者 均 在 一 定 程 度 上 参与 马尾 松 对 非 生物 胁迫 的 响应 〈 李 荣 ，2021; 朱 灵 芝 ，2021) 。 
PmHDR 蛋白 可 能 在 一 定 程 度 上 促进 马尾 松 的 生长 发 育 和 逆境 响应 过 程 。 本 文 研究 发 现 PmHDR 基因 

长 度 为 1 458 bp， 编 码 485 个 氨基 酸 ， 其 编码 蛋白 包含 LytB/IspH 基因 超 家 族 的 核心 序列 和 PLN02821 多 功 


能 结构 域 ， 属 于 1- 羟 基 -2- 甲 基 -2-(E)- 丁 烯 基 -4- 焦 磷酸 还 原 酶 家 族 ， 这 说 明 PmHDR 蛋白 参与 催化 熙 类 


化 合 


物 的 生物 合成 。 对 PmHDR 蛋白 进行 结构 预测 ， 发 现 该 蛋白 富 仿 a 螺 旋 与 B 折 辣 ， 三 级 结构 呈 “ 首 落叶 ” 形 ， 


与 其 他 物种 中 的 HDR 蛋白 结构 相似 ，4 个 保守 的 半 胱 氨 酸 残 基 可 能 处 于 “ 音 蒂 叶 ” 中 间 位 置 ， 构 成 铁 硫 建 


形成 酶 的 催化 中 心 ， 从 而 参与 维持 细胞 内 的 氧化 还 原 反 应 《〈 刘 苗 苗 等 ，2017) 。 


研究 不 同 植物 HDR 基因 表达 和 功能 时 ， 需 要 根据 物种 特性 及 其 主要 碳 类 化 合 物 的 类 型 和 功能 ， 有 针 
对 性 地 进行 分 析 和 验证 。 本 研究 结果 发 现 PmHDR 基因 在 马尾 松 各 个 组 织 中 均 有 表达 ， 但 在 老 叶 中 表达 量 


明显 高 于 其 他 部 位 ， 其 次 为 幼 叶 、 幼 茎 和 老 茎 ， 在 根 中 表达 量 最 低 ， 这 与 铁皮 石 解 ( 吴 秋 菊 等 ，2015) 和 
H3 (Dioscorea esculenta) (Wang et al., 2012) 中 HDR 的 表达 特征 较为 相似 。 然 而 在 葛 革 木 (Chen et al., 
20100 F, HDR 基因 在 花 中 表达 量 最 高 ， 果 实 次 之 ; 在 秦 欧 (Gentiana macrophylla) (XÆ, 2015) rp, 
HDR 基因 在 花 中 表达 水 平 显 著 高 于 其 他 组 织 ， 而 在 根 、 茎 、 叶 中 表达 差异 不 显著 ; 在 银杏 (u etal., 2008) 
和 橡胶 (Sando et al., 2008) 中 ，HDR 基因 在 根 中 表达 量 则 显著 高 于 茎 和 叶 。 这 说 明 HDR 基因 在 不 同 植物 
中 的 表达 模式 和 表达 量 存在 较 大 的 差异 ， 这 可 能 是 因为 在 不 同 植物 中 ， 次 生 代谢 所 产生 的 莫 类 化 合 物 类 型 
和 数量 不 同 ， 进 而 导致 其 功能 和 应 用 的 差异 ， 因 此 对 于 不 同 物种 中 HDR 等 帖 类 化 合 物 生 物 合成 上 游 关 键 


基因 的 研究 目的 也 各 不 相同 。 如 对 于 秦 苑 〈 崔 文 等 ，201$) 、 艾 叶 〈 刘 苗 苗 等 ，2017) . FF (Chen etal., 


2013) 等 药 用 植物 来 说 ， 页 类 及 机 类 衍生 物 是 其 重要 的 有 效 药 用 成 分 ， 对 于 葡萄 (Vitis vinifera) (Chen et 
al., 2023) 来 说 ， 苘 烯 类 物质 是 葡萄 香气 的 主要 贡献 者， 而 香气 是 衡量 葡萄 品质 的 重要 指标 。 

对 密码 子 偏 好 性 进行 分 析 ， 不 仅 有 助 于 基因 功能 、 和 蛋白 表达 的 研究 ， 亦 有 助 于 选择 合适 的 宿主 ， 提 高 
外 源 基 因 表 达 效 率 ( 吴 宪 明 等 ，2007) 。 本 研究 对 马尾 松 PmHDR 基因 的 密码 子 偏 好 模式 进行 了 系统 分 析 ， 
发 现 PmHDR 基因 在 编码 蛋白 时 对 CCA、AGA、GGA、ACU 和 CUU 等 5 个 密码 子 的 偏好 性 极 强 ， 而 编码 
RARAJ CCC, CCG 和 编码 Arg 的 CGG 密码 子 可 能 不 参与 PmHDR 基因 的 翻译 过 程 。ENC、CAI、FOP 
等 多 个 参数 显示 ，PmHDR 基因 对 密码 子 的 选择 偏好 性 弱 ， 这 可 能 导致 其 编码 蛋白 表达 水 平 较 低 。PmHDR 
基因 中 密码 子 的 GC 和 GC3s 含量 分 别 为 39.92% 和 34.77%， 说 明 其 密码 子 偏 好 以 AT 结尾 ， 这 一 结果 符合 
马尾 松 总 体 上 偏好 使 用 第 3 位 为 A/T 碱 基 的 密码 子 的 特征 (Zhu et al., 2020) o 
由 于 马尾 松 同 源 转化 平台 的 搭建 存在 难度 ， 至 今 未 完成 搭建 ， 目 前 的 功能 基因 分 析 只 能 暂时 借助 异 源 
转化 。 本 研究 将 PmHDR 基因 的 密码 子 使 用 频率 分 别 与 拟 南 芥 、 烟 草 、 酿 酒 酵母 和 大 肠 杆 菌 等 模式 物种 的 
基因 组 进行 比较 分 析 ， 发 现 相 比 于 大 肠 杆菌 原核 表达 系统 ， 酵 母 真 核 表 达 系 统 更 适合 作为 PnHDR 基因 的 
表达 系统 ， 而 在 PmHDR 遗传 转化 功能 验证 中 ， 拟 南 芥 和 烟草 均 可 作为 其 遗传 转化 受 体 ， 这 一 结果 与 在 
PmDXR 基因 中 的 研究 一 致 (Zhu et al., 2021) 。 因 此 在 进行 PmHDR 基因 功能 验证 时 ， 本 研究 使 用 了 拟 南 
芥 作为 其 遗传 转化 受 体 ， 通 过 对 比 转基因 拟 南 芥 和 野生 型 拟 南 芥 在 和 干旱 和 盐 胁迫 条 件 下 的 表 型 差异 ， 探 究 
该 基因 的 功能 。 
在 正常 生理 条 件 下 ， 转 基因 拟 南 芥 生长 状况 略 好 于 野生 型 拟 南 芥 ， 但 差异 并 不 显著 。 在 干旱 和 盐 胁 迫 
处 理 条 件 下 ， 转 基因 拟 南 芥 生 长 受 抑制 情况 显著 弱 于 野生 型 拟 南 芥 ， 表 现 出 更 强 的 抗 道 性 ， 这 与 其 上 游 基 
因 PmDXS (XR, 2021) 和 PmDXR CEX E, 20210 的 研究 相似 。 前 人 研究 已 经 发 现 ， 在 葛 鞭 木 (Chen 
etal,20100 和 铁皮 石 狸 〈 吴 秋 菊 等 ，2015) F, HDR 基因 的 表达 均 受 ABA 的 调控 。 作 为 一 种 胁迫 激素 ， 
ABA 在 植物 抗 逆 过 程 中 起 着 重要 人 作用， 广泛 地 参与 到 植物 抗 寒 、 抗 旱 、 耐 盐 等 多 个 胁迫 响应 过 程 。ABA 
的 生物 合成 主要 有 类 机 途径 、 类 胡 葛 小 素 途 径 两 条 途径 ， 在 高 等 植物 中 以 类 胡 葛 小 素 途 径 为 主 〈Yoshida et 
N al., 2020) , iu2S 8H b ze iX MEP 途径 的 下 游 通 路 ， 因 此 PmHDR 基因 对 干旱 和 盐 胁 人 据 的 响应 很 
c2 可 能 与 ABA 有 关 ， 且 可 能 存在 反馈 调节 。 这 一 结果 可 以 为 马尾 松 抗 旱 和 抗 盐 碱 的 分 子 育种 提供 候选 基因 
nu 和 参考 ， 但 PmHDR 基因 调控 干旱 和 盐 胁 迫 响 应 的 具体 调节 机 制 还 需要 进一步 研究 。 
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